Self-consistent calculations of the electric giant dipole resonances 
in light and heavy mass nuclei. 
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While bulk properties of stable nuclei are successfully reproduced by mean-field theories employing 
effective interactions, the dependence of the centroid energy of the electric giant dipole resonance 
on the nucleon number A is not. This problem is cured by considering many-particle correlations 
beyond mean-field theory, which we do within the Quasiparticle Time Blocking Approximation. 
The electric giant dipole resonances in ie O, 40 Ca, and 208 Pb are calculated using two new Skyrme 
interactions. 

PACS numbers: 21.30.Fe, 21.60.-n,24.30.Cz, 21.10.-k 



The electric giant dipole resonance (GDR) is a well- 
known nuclear excitation mode which is related to bulk 
properties of nuclei, such as the Thomas-Rcichc-Kuhn 
(TRK) sum rule and the nuclear symmetry energy 
One might assume that theories which describe both bulk 
properties of nuclei and shell effects rather well, such 
as self-consistent mean-field theories based on effective 
nucleon interactions 0^H, should have no problem in 
systematically reproducing the centroid energies of the 
GDR as a function of the nucleon number A. This is 
not the case, however, as has been discussed in detail 
in several recent reviews on mean-field theories which 
include strength functions obtained within the quasi- 
particle random-phase approximation (QRPA) [6|49|. It 
was impossible so far to describe ground-state properties 
and the centroid energy of the GDR both in light and 
heavy nuclei with the same effective interaction. The 
problem is more serious than might appear at a first 
glance because the physics of the GDR is intimately re- 
lated to the neutron skin thickness and the pygmy dipole 
strength T3-12|, presently investigated experimentally 
because of an impact on the isotope abundance produced 
in supernova explosions [13}. There are two hints sug- 
gesting that the mean field approach by itself is at the 
origin of the problem. Complex configurations play a 
well-known role in the damping of nuclear excitations 
14 1 . Even when effective interactions are fitted to the 
effective isoscalar mass, the symmetry energy, and the 
TRK sum rule enhancement factor k, the problem re- 
mains unsolved Q. 

We employ the Quasiparticle Time Blocki ng App rox- 
imation (QTBA), developed and applied in |15l421|. to 
study the GDR. The QTBA is a method to calculate 
nuclear response functions which generalizes the QRPA. 
It includes explicitly the coupling of one-particle one- 



hole(lplh) configuration with phonons, but omits the 
simultaneous excitation of two-phonon states in the pres- 
ence of a lplh-excitation. In the limit of vanishing 
phonon-nucleon coupling, the QTBA corresponds to the 
QRPA, a standard mean field approach. Originally, the 
QTBA was used in the framework of Landau- Migdal the- 
ory, but has been generalized recently to effective interac- 
tions of the Skyrme family in order to make possible self- 
consistent calculations O 0, 22 1- The Skyrme inter- 
actions are denned by a set of momentum- and density- 
dependent contact interactions; different parameteriza- 
tions may be distinguished by some set of theoretical 
quantities, such as nuclear matter properties or the ef- 
fective mass, which are not directly observable. The mo- 
mentum dependence of the Skyrme interaction leads to 
an effective mass, with values m*/m < 1 found by many 
investigations. Mean- field approaches which employ ef- 
fective masses smaller than unity generate single-particle 
energies which systematically deviate from the separation 
energies, mainly by a too small level density. Larger level 
densities can be obtained by taking into account the en- 
ergy dependence of the nucleon self-energy, as is shown in 
Refs. [23J426J. The energy dependence of the self-energy 
is due to complex configurations, such as the coupling 
of phonons to the single particle degrees of freedom. In 
this communication we show that if these effects are con- 
sidered, both centroid energies and spreading widths of 
the giant resonances are reproduced. As we know that 
Skyrme forces cannot reproduce simultaneously the GDR 
in 16 and 208 Pb ||, we have adjusted new Skyrme pa- 
rameterizations for the purpose of this study concentrat- 
ing on tuning the GDR in 16 O within the mean- field ap- 
proach (RPA). Since there are only few collective nuclear 
vibrations in light nuclei, the inclusion of phonons within 
the QTBA is expected to produce results close the ones 
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obtained in the mean field approach for 16 0. On the other 
hand, in heavy nuclei, the number of collective modes 
increases, which leads to major differences between the 
mean-field approach and the QTBA. We follow exactly 
the same fitting strategy and data as used for the system- 
atic variation of forces in 0. As a result we obtain two 
new forces, SV-m56-0 with effective mass m* /m = 0.56 
and SV-m64-0 with m*/m = 0.64. Both forces have a 
rather low symmetry energy Osy m = 27 MeV, and high 
sum-rule enhancement factor |6j ktrk = 0.6. The pa- 
rameters are listed in Table HI 

TABLE I. Skyrme force parameters (upper block), adjusted 
nuclear matter properties (middle block), and dipole polariz- 
ability ao as well as neutron skin r rms , n — r rmS)P (lower block) 
for the two newly designed Skyrme forces. The standard force 
parameters are given where a is the power of the density de- 
pendence. All three forces use Coulomb exchange in Slater 
approximation and the cm. energy correction (P^ vl )/(2mA). 
For details of the functional and options see [H, S0]- 





SV-m56-0 


SV-m64-0 


to 


-1905.403 


-2083.855 


tl 


571.187 


484.604 


t 2 


1594.803 


1134.345 


t 3 


8439.036 


10720.663 


U 


133.268 


113.973 


x 


0.644020 


0.619768 


Xl 


-2.973738 


-2.332678 


X2 


-1.255261 


-1.305938 


X3 


1.796625 


1.210109 


b' 4 


52.97011 


62.92567 


a 


0.2 


0.2 


ft 2 

Imp 


20.74982 


20.74982 


h 2 
2??i„ 


20.72126 


20.72126 


m* /m 


0.56 


0.64 


asym/MeV 


27 


27 


KTRK 


0.6 


0.6 


«d /fm' 1 


20.2 


19.4 


n-skin [fm] 


0.156 


0.134 



Of course, the new fits maintain the good ground state 
properties of all the systematically varied forces in [6]. 
Additionally, the low effective mass and low symmetry 
energy a sym =27 MeV together with a rather high sum- 
rule enhancement factor ktrk of the two parameteriza- 
tions delivers a high GDR energy. This is beneficial for 
16 but leads to its overestimate in 208 Pb. However the 
inclusion of complex configurations brings the GDR in 
208 Pb down to the correct value. Similar effects are seen 
for the single-particle energies. The parameterizations 
SV-m56-0 and SV-m64-0 reproduce these energies in 
16 O reasonably well, but the single-particle spectrum in 
208 Pb is spread out too much and deviates strongly from 
the experimental one. The coupling to the phonons will 
improve the spectrum Being built on RPA, QTBA 
follows basically the same trends with varying Skyrme 
force as RPA (see [6]), though the absolute value of the 



effect is different in light and heavy mass nuclei. 

Let us outline some technical details of our numerical 
scheme. In our RPA and QTBA calculations of the GDR, 
the single-particle continuum is treated exactly according 
to the scheme described in Ref. The phonons were 

calculated within the so-called discretized RPA (DRPA). 
Here the results depend on the single particle basis and on 
details of the discretization e.g. the size of the box one 
chooses. In the present investigation, such ambiguities 
are small, as we control them by comparing the DRPA 
results with a full continuum RPA. 

In self-consistent calculations the p/i-intcraction is 
given by the second derivative of the energy functional. 
In the DRPA calculations of phonons the matrix ele- 
ments of the ph-interaction were calculated exactly ex- 
cept for the spin-orbit and Coulomb contributions which 
were omitted. In the RPA and QTBA calculations of the 
GDR we used additional (local-exchange) approximation 
for the velocity-dependent part of the ph-interaction de- 
rived from Skyrme energy functional. In the case of the 
GDR this approximation gives results which are close to 
the exact RPA results. It will be described in a forthcom- 
ing publication. In our investigation we did not study 
spin-dependent properties of new Skyrme forces. It is 
well-known that spin-spin part of the residual interaction 
(except for J 2 -generating terms) does not contribute to 
the ground-state structure for spherical even-even nuclei. 
So one can omit this part in the calculations of the excited 
states in these nuclei without breaking self-consistency. 
On the other hand, inclusion of the spin-spin part of 
the residual p/i-interaction leads to an instability in the 
DRPA calculations of the phonon's characteristics with 
given Skyrme forces. In fact, as for most other Skyrme 
forces, the instability is driven by the term cx (Vcr) 2 in 
the functional. For this reason we exclude this part of 
the interaction in our calculations. 

The number of solutions of the RPA equations depends 
on the size of the configuration space. However, the 
majority of the RPA wave functions are dominated by 
one particle-hole configuration. In principle, one has to 
subtract the second order contributions to complex (ph- 
phonon) configurations in order to avoid double counting 
[27J . For simplicity, the present calculations consider only 
a small number of phonons, defined by having transition 
probabilities of at least 1/5 of the strongest state of each 
multipolarity. For these phonons, the second-order cor- 
rections are small and have been neglected. 

In FigJTJ we show the sensitivity of the photo absorp- 
tion cross sections, obtained in the framework of RPA, 
on small variations of the effective mass. We used two 
different values m*/m=0.56 and 0.64. The higher ef- 
fective mass gives lower GDR energies in all three nu- 
clei. As we are interested in a Skyrme paramctrization 
which reproduces the GDR in 16 O we present here only 
results for the lower effective mass. Note that the present 
calculations improve the description of the GDR in 16 O 
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FIG. 1. (Color online) Photo absorption cross-section 
in ie O calculated self-consistently in RPA, using two dif- 
ferent Skyrme parameterizations with effective mass 0.56 
(dashed(blue)) line) and 0.64 (dashed-dotted(green) line). 
The experimental cross section are given by the (brown) dots 
connected by a solid line [28l |. 



FIG. 2. (Color online) Comparison of the experimental [28|] 
photo absorption cross-section in 16 O with theoretical ones 
calculated in RPA (dashed(blue) line) and QTBA (solid(red) 
line) using the SV-m56-0 Skyrme parameters. The experi- 
mental data are given by the (brown) dots connected by a 
solid line 



TABLE II. Comparison of theoretical and experimental [2^ - 
[30| Lorentzian parameters. The energies considered range 
from 8-25 [MeV] for 208 Pb and from 10-32 [MeV] for 40 Ca 
and 16 O. 



Nucleus 


Force 


E [MeV] 


r [MeV] 


(Jo [mb] 




SV-m56-0 (RPA) 


14.30 


4.96 


624 




SV-m56-0 (QTBA) 


13.37 


5.99 


495 




Experiment 


13.43 


5.08 


481 


4l) Ca 


SV-m56-0 (RPA) 


21.61 


5.90 


104 




SV-m56-0 (QTBA) 


21.14 


5.92 


99 




Experiment 


20.00 


5.00 


95 




SV-m56-0 (RPA) 


25.31 


8.95 


25.7 




SV-m56-0 (QTBA) 


24.49 


8.85 


24.8 




Experiment 


23.76 


7.17 


24.8 



800 




E [MeV] 



in comparison with other self-consistent approaches. As 
only a few collective states exists in 16 O, we do not ex- 
pect strong modifications of the RPA results due to the 
phonons. This is indeed the case and is demonstrated 
in Fig. [2] where we compare the RPA and QTBA results. 
The Lorentzian parameters of the photo absor ptio n cross- 
section [l[ derived from the data of Refs. |28l43cl ] are 
summarized in Table E] for 208 Pb, 40 Ca, and 16 0. The 
data shown here and in the subsequent figures [28l - [30l | 
are also available electronically (3lj. In Fig. [3] we present 
the dipole photo absorption cross-section in 208 Pb calcu- 
lated with the Skyrme parametrization SV-m56-0 and 
compare them with the data [29| . The result of the con- 
ventional RPA is compared with the QTBA where the 
phonons are included. In RPA the mean energy of the 
GDR E= 14.30 MeV is too high. The rather large width 
r=4.96 MeV in the RPA is explained by a strong peak 



FIG. 3. (Color online) Comparison of the experimental photo 
absorption cross-section [2ij] in 208 Pb ((brown) dots with 
bars), with theoretical ones calculated in RPA (dashed(blue) 
line) and QTBA (solid(red) line). The Skyrme parametriza- 
tion SV-m56-0 was used. 

in the cross-section at 20.4 MeV, i.e. in the high-energy 
tail of the GDR. The phonons shift the GDR to lower en- 
ergies, where the mean energy E = 13.37 MeV and the 
width T=5.99 MeV are now in good agreement with the 
experimental data. We investigated the photo absorp- 
tion cross section in 40 Ca with the same Skyrme force 
SV-m56-0 as an example for an intermediate mass nu- 
cleus which is shown in Fig. |4] The RPA result is about 
1.6 MeV higher compared to the data. The cross-section 
calculated within the QTBA is shifted by 0.5 MeV to 
lower energies and agrees better with experiment. 
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FIG. 4. (Color online) The same as in Fig. 02 but for 40 Ca. 
The data are taken from Ref. 13011. 



In summary, we show that the explicit inclusion of 
quasiparticlc-phonon coupling may solve the problem of 
mean-field theories in reproducing the centroid energies 
of the giant dipole resonance. As the phonon contribu- 
tion is small in light nuclei, but large in heavy mass nu- 
clei, phonon excitations provide a mass-dependent mech- 
anism for damping and energy shift. First calculations 
employing two new Skyrme interactions show a reason- 
able quantitative agreement with the experimental dipole 
excitations in 16 , 40 Ca, and 208 Pb. 
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